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ABSTRACT 

We study higher order corrections to Higgs production with an associated jet 
at SSC energies, using the resummation of the leading logarithmic contributions to 
multiple gluon emissions due to Lipatov and collaborators. We find a considerable 
enhancement of Higgs production at large transverse momenta. 
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1. Introduction 



One of the most important open questions in the Standard Model is the cause 
of the electroweak symmetry breaking. Both the SSC and the LHC are being 
buih with the resolution of this puzzle in mind. The simplest explanation of the 
symmetry breaking is by a fundamental SU{2) doublet scalar field, which leaves 
a neutral scalar, the Higgs boson, as a remnant of the symmetry breaking. This 
simple theory is the experimental standard against which any model of symmetry 
breaking will be measured. Thus, it is important that we completely understand 
the production environment of the Higgs boson at these hadron colliders. 

At the SSC and the LHC, the Higgs boson is predominantly produced by 
gluon fusion. At lowest order a^, it is produced with zero transverse momentum. 

o [11 

The complete order corrections have been calculated with an increase in the 

total cross section by a factor of 1.5 to 1.7. In addition, production of Higgs 

bosons with nonzero transverse momentum have been studied at order cc^. Higgs 

bosons produced with sizable m^-y be phenomenologically important both in the 

intermediate mass region, where rare decay modes such as 77 and t'^t~ must be 

used to observe the Higgs, and for heavier Higgs bosons where the vector boson 

[3] 

decay channels may be used. At small transverse momentum it is necessary to 
resum large logarithms in m'jj/p'j_ in order to obtain a physical p± distribution.^^ 

At the large center of mass energy y/s of the SSC and the LHC there are also 
important contributions from events with multiple final-state partons. These arise 
in a new kinematical region, the semihard region, characterized by scattering pro- 
cesses with y/s much larger than the momentum transfer Q, s >> >> ^qcd- 
The study of this region is theoretically challenging because, in the perturbative 
expansion of cross sections, there appear coefficients containing logarithms of large 
ratios of kinematical invariants, of the order of the rapidity interval in the scat- 
tering process. In the inclusive production of the Higgs boson, the large ratios of 

kinematical invariants appear in the small-x evolution of the structure functions, 

[7] 

thus requiring a more sophisticated analysis than the usual DGLAP evolution. 
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In the context of jet production, Mueller and Navelet were able to disentangle 
the small- a; evolution of the structure functions from the appearance of large ratios 
of kinematical invariants in the parton subprocess by proposing to tag two jets at 
the extremes of the Lego plot in azimuthal angle and rapidity, with large transverse 
momenta pj^, and longitudinal momentum fractions xa, xb large enough that the 
parton structure functions could be computed from the ordinary DGLAP evolution. 
Thus, fixing the parton center of mass energy § = xaxbs, the semihard region is 
realized also at the parton level s >> Q^, and the large logarithmic terms \n{s/Q'^) 

appear in the parton cross section. To deal with these, Mueller and Navelet used 

[91 

the Balitsky-Fadin-Kuraev-Lipatov theory (BFKL), which systematically resums 
the leading powers in the rapidity interval by using a multigluon amplitude where 
the rapidity interval between the tagging jets is filled with gluons, the minijets, 
whose spacing in rapidity is approximately uniform. 

For the Higgs boson this program can be implemented by considering the pro- 
duction of the Higgs boson and a jet at the extremes of the Lego plot, with the 
rapidity interval between them filled with minijets, and the Higgs boson and the 
jet tagged at large transverse momenta. Since in gluon fusion the Higgs boson is 
produced via a quark loop, and the fermion-Higgs coupling is proportional to the 
quark mass, the only contribution to the quark loop which is numerically impor- 
tant is due to the top quark. Thus, two more scales, the Higgs boson and the 

top quark masses, enter the kinematics, and in order to avoid complications with 

[7] 

the evolution of the structure functions and apply the BFKL analysis, we must 
require that s >> |)^,m^,m^. However, when « doubly logarithmic 
terms ln^(m|^/p^) appear, which more properly belong to the evolution of the 
structure functions and may void the BFKL analysis. So in order to single out 
the strong rapidity-ordering regime characteristic of the BFKL analysis, the Higgs 
and the jet transverse momenta must be comparable to the Higgs mass ~ mjj. 

In section 2, we consider the inclusive production of the Higgs boson and a 
tagging jet, pp ^ H + jet + X in this kinematic situation. We recall the exact 
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Born-level calculation of Refs. 2 and 3 and consider its large-rapidity limit, and 
then we compute the leading logarithmic corrections at large rapidity, using the 
BFKL analysis. In section 3, we present numerical results for the inclusive Higgs-jet 
production, with or without integrating over the Higgs and jet transverse momenta, 
and compare the Born-level to the minijet-corrected calculations. Our conclusions 
are presented in section 4. 

2. The Higgs + Jet Inclusive Cross section 

We are going to study the semi-inclusive process pp ^ H + jet + X in the 
semihard regime defined by s >> Q^, with being a typical momentum scale in 
the event, ~ m'jj,m^,p'^. The tagging jet is required to ensure that we have 
an event with a large rapidity interval y = yj — yu ~ ln(s/Q^). Other relevant 
parameters in the event are the Higgs and jet transverse momenta Ph±,PJ±, their 
relative azimuthal angle 4> and the average rapidity y = {yj + yH)/2. 

In the semihard, large-y regime we can write the cross section: 



The parton subprocess cross section daij/dp'jjj_dpjj_d(f) contains the sum over all 
additional particles (i.e. minijets) in the event. The factorization of the minijets 
into the subprocess cross section is possible, because at large y the initial parton 
momentum fractions xi and X2 are fixed in terms of the Higgs and jet momenta, 
and are essentially independent of the particles filling the rapidity interval. We will 
arrive at this cross section in several steps, starting with the exact Born level cross 
section, taking it to the y » 1 limit, and finally filling in the rapidity interval 
with the minijets. Of course, the analysis can be applied equally well for y << — 1, 
by reflecting along the beam axis. 

i) Born Level Cross section . At the Born level the Higgs boson and the parton 
jet are produced back-to-back. The exact lowest order cross section can be put in 




(2.1) 
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the form (2.1) with the replacement 



da. 



dt 



(2.2) 



The parton momentum fractions and the subprocess invariants at this level are 
given by: 

ey 



X2 = 

s — 

i = 

u — 



XIX2S 

-pI- 

-pI- 



p\ + m]j^ + 2p_LmH± cosh(y) 
,-y 



(2.3) 



P±rnH±e 



P±mH±e 



y 



where mui. — {p\ + m?jj)^f'^ and p±^ — pui. — Pja_- We use the lowest order 
Standard Model calculation for the subprocess cross sections daij/di in Ref. 2. 

■i/i) Largc-y Born Cross section. We now investigate the lowest order cross sec- 
tion when the rapidity interval y is large. For y » 1 the lowest order amplitude 
is dominated by diagrams with gluon-exchange in the t-channel as in Fig. 1(a). In 
this limit the only subprocesses that contribute are gg — > gH and q{q)g q{q)H. 
We obtain 



da, 



gg 



dt 



C 



\Hp 



pI 



where 



C = 



Ncodaw 



(7V2 _ i)i287rM2, " 



Similarly, we find 



da, 



qg 



da-qg 



Cp da. 



99 



(2.4) 



(2.5) 



(2.6) 



dt dt Ca dt 

with Cf/Ca = {NI — 1)/2NI = 4/9 the ratio of the Casimir operators. The form 
factor J-{p\) is given in the appendix for both a scalar and a pseudoscalar Higgs 
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boson. We have checked that (2.4) agrees with previous resuhs if we take the 
additional Umit p\ « mjj,m^. The parton momentum fractions in the large-y 
hmit are 

Equation (2.7) is also valid in the large- y limit when higher-order corrections are 
included, so that pui. 7^ 

ill) Minijet-corrected Cross section . As discussed in the introduction, going 
to higher orders in the coupling constant, i.e. to multiple parton emission, we 
encounter large logarithmic contributions. In the semihard regime, the BFKL 
theory systematically resums the leading logarithmic terms ln(s/p^) by using a 
multigluon amplitude where the rapidity interval between the Higgs boson and the 
tagging jet is filled with gluons, strongly ordered in rapidity. This amplitude is 
shown in Fig. 1(b), where the thick line represents the resummation of the virtual 
radiative corrections, whose effect is to reggeize the gluons exchanged in the t 
channel. The real gluons are inserted on these using the Lipatov effective 3-gluon 

[9] 

vertex. The BFKL multigluon amplitude is then put in a rapidity-ordered phase 
space, the rapidities of the gluons are integrated out, and the dependence of the 
cross section on the gluon transverse momenta is reduced to the resolution of an 
integral equation. Its solution is then convoluted with the Higgs-boson production 
vertex on one side and the jet emission vertex on the other side of the rapidity 
interval to give the minijet-corrected parton cross section for producing the Higgs 
boson and a jet: 

da _ 2C\T{pIJ\\ 

j[y,PH±,pj±) ■ (2.sj 



In this equation f{y,PH±,Pj±) is the Laplace transform in the rapidity interval y, 

^^e'^^L{PH±,PJl.), (2.9) 
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of the solution of the BFKL integral equation 



OO 

1 r l'r)2 \-l/2+ii'(J2 \-\l2-w 

Uph^,Pj^) - j:^, E ^^'^^'-'^^ / d^^^ H ■ (2.10) 

ZTT r ^-^ J u} — u}{n,u) 

^ ' n=—oo \ ' / 



The eigenvalue of the integral equation a;(n, u) is 



u;(n, v) = [^(1) - i?e V'l^ri^ + iv)] , (2.11) 

with %l) the logarithmic derivative of the Gamma function. Substituting (2.9) and 
(2.10) in (2.8), and doing the integral over u, the minijet-corrected parton cross 
section becomes 

OO 

da ^ 1-^^-2 M2 . / ^2 



^^l^PHjyLy e^ni^-) f di^e^^-'^^y COS (i^ln^A 



(2.12) 

If we integrate over the azimuthal angle in (2.12), only the n = term survives. 

iv) Minijet-corrected Cross section in the Saddle Point Approximation. 
At very large values of the rapidity interval y, the correlations between the Higgs 
boson and the jet are washed out by the random walk in transverse momentum 
space of the gluons exchanged in the t channel. This can be seen most easily by 
evaluating (2.12) in the saddle-point approximation. The contribution of (2.11) to 
this equation is dominated by n = and is strongly peaked near v — ^. Thus we 
keep only the first term in the Fourier expansion in 0, and expand Lo{y) = a;(0, v) 
about i/ = 



with 



= + (2.13) 



— ^ln2, S = ^C(3). (2.14) 

TT TT 



Then we can evaluate (2.12) using the saddle-point approximation for the integral 
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over z/, to obtain 



da C e^y 



exp 



( 



ABy j PH±P^ji_ 



(2.15) 



The exponential growth of (2.15) with the rapidity interval y is due to the produc- 
tion of the minijets. 



We now examine numerically the effects of the minijets on Standard Model 
Higgs production at SSC energies. Throughout this section we will use a repre- 
sentative Higgs mass of 100 GeV and a top mass of 150 GeV, and the SSC center 
of mass energy ^/s — 40 TeV. The qualitative results do not depend strongly on 
these values of mu and m^. Except where indicated we set y = and observe 
the cross sections as a function of the rapidity interval \y\. We obtain a factor 
of two by including both positive and negative y. We use the "average" parton 
density functions of Diemoz et al.]^^^ and we evaluate the QCD coupling ag and 
the structure functions at a scale = mj^. For pj_ > 50 GeV and > 4 the 
structure functions are always evaluated at x > 10~^, so we are justified in using 
the DGLAP evolution in this region of phase space. 

In Fig. 2 we present the inclusive Higgs- jet cross section, integrated over the 
azimuthal angle (f) and over both transverse momenta with a cutoff of p_Lmin = 50 
GeV. We present it in the four approximations (i-iv) given in section 2. From the 
plot we see that the large-y Born cross section is a good approximation to the exact 
Born level cross section for \y\ ^ 4. The minijet-corrected cross section is typically 
enhanced over these lowest order cross sections by a factor of 2-3 for \y\ ^ 4. We 
define the i^-factor as the ratio of the minijet-corrected cross section to the large-y 
Born cross section 



3. Numerical Results 




(3.1) 



and plot it in Fig. 3. K is defined so that K 



1 as y 



0. The minijet 
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enhancement is often plotted on a logarithmic scale ; on such a plot, the exact 
integral approaches the saddle point approximation, since the relative error tends 
to zero at large s. 

We next show the p± distributions of the Higgs and the tagging jet for rapidity 
intervals of \y\ = 4 and \y\ = 8 in Figs. 4(a) and 4(b), respectively. The other 
momentum is integrated over pj_ > p±rnm = 50 GeV. Of course, for the Born cross 
section we have pni. = Pj±, so the transverse momentum distributions of the 
Higgs and the jet are identical. However, with the inclusion of the minijets, this is 
no longer true. In fact the pjji. distribution is considerably flatter than either the 
lowest order result or the minijet-corrected pj± distribution. 

This effect can be understood by referring to Fig. 1(b) and analyzing the var- 
ious terms in the minijet-corrected cross section (2.8). When convoluted with the 

structure functions we find 



fiy,PH±,Pj±) 



/ 1 \2' 

xif{xi)[ ) 



(3.2) 



where the first bracketed term is associated with the Higgs boson production vertex 
and the last bracketed term is associated with the jet emission vertex. Note that 
the comes from the jet vertex, whereas the form factor ^{p'h±) is relatively 

constant in the region of interest. For very small rapidities we approach the Born 
cross section (2.4) with 

f(y,PH±,PJl.) S{pjj^-p'j_^)5{(l)-7r) , (3.3) 

so that both the Higgs boson and jet pj_ distributions fall as In addition, 

both xi and X2 increase with p^ via (2.7) so that there is an additional suppression 
from both parton density factors {xf{x))'^ as increases. However, for very large 
rapidities we have 

f{y,PH±,pj±) ~ {pH±pj±y^ , (3.4) 

and the Higgs boson and the tagging jet become uncorrelated. Now the Higgs 
boson distribution falls slower as 1/phl-i while the jet distribution falls faster as 
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In addition both distributions are now only suppressed by one factor of 
xf{x). These effects combine to strongly broaden the Higgs distribution. For 
the tagging jet the two effects tend to cancel and its p± distribution is much less 
modified. 

Comparing Figs. 4(a) and (b), we can also observe some subtleties that have 
been neglected in the arguments of the last paragraph. For example, the separation 
between the Higgs and jet pj_ distributions at large pj_ shrinks as \y\ goes from 4 
to 8. This is due to the fact that both quarks and ghions can enter the parton 
scattering on the jet side, while only gluons can enter on the Higgs side. Thus, the 
suppression from the structure functions is more severe on the Higgs side than on 
the jet side as x grows according to (2.7). 

Finally, we discuss the effects of varying the Higgs and top-quark masses in the 
p± distributions of Fig. 4. Although the change in the contribution of the minijets 
is minor, there will still be some overall effects on the distributions. For example, 
because we have chosen the factorization and renormalization scales to be fixed at 
the Higgs mass, as uih grows the coupling constant ctg and the structure functions 
become smaller, while the parton momentum fraction x increases. Consequently, 
the absolute scale of the p± distributions will decrease. Changing the top-quark 
mass has a different effect. The form factor ^(pj^j) is roughly constant for pj_ < 
rrit, but falls off as m^/p^ times logarithms for p± ^ mt- Thus, while the jet p± 
distribution becomes independent of mt at very large y, the Higgs and Born p±^ 
distributions fall more steeply for pj_ ^ m,t. For a given fixed p±_ larger than the 
top mass, the Higgs and Born pj_ distributions will increase with m^. 

The disappearence of correlations as \y\ increases can be seen dramatically in 
Fig. 5 where we plot the transverse momentum distribution of the tagging jet at 
a fixed Higgs transverse momentum of ph± = 125 GeV. For a rapidity of \y\ = 2 
the cross section is strongly peaked near pij± = pjj_. As the rapidity is increased 
there is a diffusion of the jet momentum away from the Higgs momentum until the 
peak disappears completely for \y\ — 8. Note that the saddle point approximation 



10 



removes the correlation for all values of \y\. There is a similar diffusion of the 
azimuthal angle away from = tt as the rapidity increases. 

In Fig. 6 we consider the pui. distribution of centrally produced Higgs bosons. 
We fix yH — while integrating over the rapidity and transverse momentum of the 
tagging jet. The dashed hue is the lowest order Born calculation with no cuts on 
the jet rapidity, while the dotdash line has a cut of > 4. The jet transverse 
momentum is fixed at = pu^ in these Born-level curves. The solid line is the 
minijet-corrected cross section with cuts of > 4 and > 50 GeV. We apply 
these cuts in order to be confident that we are in a kinematical region where the 
minijet-corrected cross section is a good approximation. As in the previous plots 
we see that the minijet-corrected pui. distribution is much broader than the lowest 
order distribution. At pui. — 400 GeV the minijet cross section with the kinematic 
cuts is over two orders of magnitude larger than the Born cross section with cuts. 
Moreover, it is almost within a factor of 2 of the lowest order Born cross section 
without any cuts at all. 

4. Discussion and Conclusions 

In this paper we have calculated the contribution of minijets to Higgs pro- 
duction in the semihard regime. The major phenomenological consequence of the 
minijets at SSC energies is to produce a substantial enhancement of Higgs bosons 
at large transverse momenta. Although our results are strictly valid for the semi- 
inclusive process pp ^ H -\- ]ei -\- X at large y, we have seen that these events may 
still make a sizable contribution to the inclusive process pp ^ H -\-X ior large val- 
ues of phj^- The enhancement occurs because in the minijet-corrected cross section 
the falloff with pui. is much slower than the falloff with pjj,. This suggests that 
there may be a sizable contribution from events with a Higgs boson produced at 
large transverse momentum, balanced by the collective p}_ of several less energetic 
jets. These events should appear first at order in the gg — > Hgg and gq — > Hgq 
cross sections. The former has been calculated at the Born level in the large mt 
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limit by Dawson and Kauffman, who plotted the cross section with a single 
cutoff on both the outgoing gluons and the Higgs boson. It would be interesting to 
study this gg — > Hgg cross section with the Higgs and gluon transverse momenta 
varied independently in order to check for an enhancement of events at large ph±- 
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APPENDIX 



The analysis in this paper is valid both for a scalar Higgs boson and for a 
pseudoscalar Higgs boson, as would occur in two-Higgs-doublet models. The only 
change in our results is in the form factor J^(p^). We present the analytic form of 
J^{p\) in both cases below. 

i) Scalar Higgs. The scalar Higgs coupling to the top quark is —S{mt/v)Hiljilj, 
where is a normalization constant which equals 1 for the Standard Model Higgs. 
The form factor is 



s 



Ami \ 
mjjj 



+ 



m 



VbW{b) - ^W{a) 
W{bf - W{af 



where a — 1 + 4m|/pj^, b — 1 — Aml/rrijj, and we take the root \/b 
6 < 0. We have also defined the function 



W{c) 



In 
In 



2i arcsin (1/ \/l — c) , 



c < 
< c< 1 
c > 1 



(Al) 



for 



(A2) 



Note that Ts{p\) is proportional to the form factor of Ref. 2 in the appropriate 
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kinematic channel. In the hmit pj^, mn « mt, the form factor goes to the value 
J^s = -(4/3)5. 

ii) Pseudoscalar Higgs . The pseudoscalar Higgs couphng to the top quark is 
—P{'mt/v)Ail)i^5ip. In this case the form factor is 



W{bf - W{af 



(A.3) 



In the limit pj_, rnu << mt, we find J^p — 2P. 
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FIGURE CAPTIONS 



1) Higgs + jet production amplitude in the \aige-y limit at (a) the Born level 
and (b) with minijet corrections. 

2) Inclusive Higgs-jet production at the SSC, as a function of the rapidity in- 
terval y. The dashed and dotdashed hues are respectively the exact and 
large-y Born cross sections, the solid hne is the minijet-corrected cross sec- 
tion, and the dotted line is the saddle point approximation to the solid line. 
The kinematic parameters are described in the text. 

3) X-f actor, as a function of the rapidity interval y. The solid line represents the 
ratio of the minijet-corrected cross section to the large-y Born cross section, 
and the dotted line is its saddle point approximation. 

4) Higgs and jet pj_ distributions at (a) \y\ = 4 and (b) \y\ = 8. The dashed 
and dotdashed lines are respectively the pj_ distributions for the exact and 
the large-y Born cross section, for which pui. = Pja_- The solid lines are the 
jet and Higgs pj_ distributions for the minijet-corrected cross section. Notice 
that in (b) the dashed and dotdashed lines completely overlap. 

5) Jet p_L distribution at a fixed Higgs transverse momentum of 125 GeV. Prom 
top to bottom, the solid lines are the jet p± distributions for the minijet- 
corrected cross section at \y\ = 2, 4, 6 and 8. The dotted lines are the saddle 
point approximation to the solid lines. 

6) Higgs p± distribution ai yn — 0- The dashed and dotdashed hues are the 
distributions for the exact Born cross section with, respectively, no cuts on 
the jet rapidity \yj\, and \yj\ > 4. The solid line is the distribution for the 
minijet-corrected cross section with \yj\ > 4 and pj^ > 50 GeV. 
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